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(Dated: Dec 5, 2006) 

We search for the technicolor process pp — > pr/uJT — > Wtvt in events containing one electron 
and two jets, in data corresponding to an integrated luminosity of 390 pb _1 , recorded by the DO 
experiment at the Fermilab Tevatron. Technicolor predicts that technipions, ttt, decay dominantly 
into bb, be, or be, depending on their charge. In these events b and c quarks are identified by 
their secondary decay vertices within jets. Two analysis methods based on topological variables are 
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presented. Since no excess above the standard model prediction was found, the result is presented 
as an exclusion in the 7tt vs. pr mass plane for a given set of model parameters. 



PACS numbers: 12.60.Nz, 13.85.Rm 

Technicolor (TC), first formulated by Weinberg and 
Susskind [l], 0], provides a dynamical explanation of 
electroweak symmetry breaking through a new strong 
SU(Ntc) gauge interaction acting on new fermions, 
called "technifermions." Technicolor is a non-Abelian 
gauge theory modeled after Quantum Chromodynamics 
(QCD). In its low-energy limit, a spontaneous breaking 
of the global chiral symmetry in the tcchnifcrmion sector 
leads to electroweak symmetry breaking. The Nambu- 
Goldstone bosons produced in this process are called 
technipions, ttt, in analogy with the pions of QCD. Three 
of these technipions become the longitudinal components 
of the W and Z bosons, making them massive. 

An additional gauge interaction, called extended tech- 
nicolor , couples standard model fermions and tech- 
nifermions to provide a mechanism for generating quark 
and lepton masses. By limiting the running of the tech- 
nicolor coupling constant, walking technicolor [H[ avoids 
flavor-changing neutral currents. To generate masses 
as large as the top quark mass, another interaction, 
topcolor, seems to be necessary, thereby giving rise to 
topcolor-assisted technicolor models [6] . 

Extensions of the basic technicolor model tend to re- 
quire the number No of teclmifermion doublets to be 
large. In general, the technicolor scale Atc w O(l) x 
Ftc, where Ftc is the technipion decay constant, de- 
pends inversely on the number of technifermion dou- 
blets: F TC ~ 246 GeV '/y/NE- For large N D , the low- 
est lying technihadrons have masses on the order of few 
hundred GeV. This scenario is referred to as low-scale 
technicolor [7|. Low-scale technicolor models predict the 
existence of scalar technimesons, iri^ and w T , and vec- 
tor technimesons, pt and lot- General features of low- 
scale technicolor have been summarized in the technicolor 
strawman model (TCSM) [g, M. The analysis presented 
in this paper is based on Ref. 0. 

Vector technimesons are expected to be produced with 
substantial rates at the Fermilab Tevatron Collider via 
the Drell-Yan-like electroweak process pp — > pt + X or 
lot + X . In walking technicolor, it is expected that vec- 
tor technimesons decay to a gauge boson (7, W, Z) and 
a technipion or to fermion-antifermion pairs. The pro- 
duction cross sections and branching fractions depend 
on the masses of the vector technimesons, M{pt) and 
M{lot), on the technicolor-charges of the technifermions, 
on the mass differences between the vector and scalar 
technimesons, which determine the spectrum of accessi- 
ble decay channels, and on two mass parameters, Ma for 
axial- vector and My for vector couplings. The parame- 
ter My controls the rate for the decay pt,lot — ► 7 + tit 
and is unknown a priori. Scaling from the QCD decay 



p, lo — > 7 + 7T°, the authors of Ref. [9( suggest a value of 
several hundred GeV. We set Ma = My, and evaluate 
the production and decay rates at two different values: 
100 and 500 GeV. For all other parameters, we use the 
default values quoted in Table III of Ref. [§| . Technipion 
coupling to the standard model particles is proportional 
to their masses, therefore technipions in the mass range 
considered here predominantly decay into bb, be, or be, 
depending on their charge. 

In this Letter, we describe a search for the decay of 
vector technimesons to Wttt, followed by the decays 
W — > ev and ttt — > bb, be, or cb. In the DO detector, 
which is described in detail in Ref. [l(J, the signature of 
this process is an isolated electron and missing transverse 
momentum (j/>t) from the undetected neutrino from the 
decay of the W boson, and two jets of hadrons coming 
from the fragmentation of the quarks from the decay of 
the technipion. Jets are reconstructed using the Run II 
cone algorithm [11] with a cone size of 0.5. We search for 
events with this signature in the data collected with a 
single electron trigger until July 2004 and corresponding 
to an integrated luminosity of 388±25 pb _1 fl2| . 

There are a number of standard model processes that 
can result in the same final state signature as Wttt pro- 
duction. Vector boson production in association with jets 
is the dominant background. Z boson production can be 
suppressed by vetoing on a second electron and requiring 
significant j>T- Most of the jets in VK+jets events orig- 
inate from the fragmentation of light quarks or gluons 
and therefore requiring the explicit identification of at 
least one jet from the fragmentation of a & or c quark 
suppresses most of this background, leaving only W + 66, 
W + b, W + cc, and W + c events. Top quark produc- 
tion followed by the decay to eub is another background. 
Top-antitop quark pair production typically results in ei- 
ther an additional lepton or a higher jet multiplicity from 
the decay of the second top quark, and this background 
can be reduced by selecting events with exactly two jets. 
Single top quark production is an irreducible background, 
but it has a smaller cross section. We simulate all these 
processes using either pythia [l]| or ALPGEN Q Monte 
Carlo (MC) generators, followed by the DO detector sim- 
ulation based on GEANT [lij]. Quark hadronization and 
fragmentation is simulated using pythia. 

The multijet background is due to events with poorly 
measured jets, resulting in missing momentum and a jet 
that is misidentified as an electron. Background from 
the mistagged IV+jets process originates from events in 
which a light-quark or gluon jet is incorrectly identified 
as a 6 jet. These instrumental background contributions 
are estimated from the same data sample before requiring 
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the identification of a b jet. 

We select events in which there is exactly one well- 
identified electron based on tracking and calorimeter data 
with transverse momentum px > 20 GeV and pseudo- 
rapidity 1 77 1 < 1.1 [16]. There must be significant ^t, 
measured in two ways: $ T 1 > 20 GeV computed as the 
negative sum of the jet momentum vectors and the elec- 
tron momentum vector and j>x > 20 GeV which also 
includes the calorimeter energy deposit not assigned to 
the electron or the jets. We require the transverse mass 
Mx{ev) > 30 GeV. We further require the presence of 
exactly two jets with px > 20 GeV and < 2.5. 

To further reduce backgrounds, we take advantage of 
the long lifetime of b flavored hadrons. Tracks from the 
decay products of b hadrons may not project back to the 
proton-antiproton collision, but have a significant impact 
parameter. They can therefore be identified and used to 
reconstruct the decay vertex of the b hadron. A jet is 
tagged as a b jet if there is a secondary decay vertex 
within ATZ = y/{A-q) 2 + (A(j)) 2 < 0.5 of the jet axis. We 
require at least one of the jets to be 6-tagged. This leaves 
us with 117 events in our final data sample. 

The expected background event yields are listed in Ta- 
ble |U When estimating these yields, each Monte Carlo 
event is weighted by the probability that at least one jet 
is tagged as a b jet. The tagging probability is param- 
eterized as a function of jet flavor, jet pt, and 77. The 
efficiency of tagging a jet from the fragmentation of a b 
quark is derived from collider data which were enriched 
in their b jet contents by requiring a muon to be recon- 
structed within at least one jet to preferentially select jets 
with semileptonic b decays. The probability of tagging 
a c jet is derived from the tagging probability for b jets 
by multiplying by the ratio of tagging probabilities for c 
and b jets derived from MC simulations. We derive the 
probability to tag a light-quark or gluon jet from a set of 
dijet events, corrected for contamination by c and b jets. 
The Monte Carlo events are also weighted by the ratios 
of jet and electron finding efficiencies in Monte Carlo and 
collider data. Electron finding efficiencies are measured 
in Z — > ee events in both data and Monte Carlo. 

We use the PYTHIA event generator to simulate sig- 
nal events, modeling initial state and final state radi- 
ation, fragmentation, and hadronization. To generate 
Witt signal events for a range of values of the technime- 
son masses, we use a fast, parameterized detector simu- 
lation that was tuned to reproduce the kinematic distri- 
butions and acceptances from events simulated with the 
detailed GEANT-based detector simulation. For the cross 
section calculations, CTEQ5L 17] parton distribution 
functions are used. Finally, as is appropriate for this 
Drell-Yan-like process, the cross section is multiplied by 
a iV factor of 1.3 to approximate NLO contributions to 
the cross section [3] ■ We generate events with pt masses 
from 160 GeV to 220 GeV and assume M(ui T ) = M(p T ). 
The 7Tt mass values start at the kinematic threshold for 



TABLE I: Number of events observed in the data and ex- 
pected from signal and background sources after the kine- 
matic selection; only statistical errors are reported. For the 
expected number of signal events quoted we assume M(pr) = 
210 GeV and M(n T ) = 110 GeV. 
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Witt production at M(7Tt) = M(px) — M(W) and go 
down to M(itt) = M(pr)/2 — 5 GeV where the decay 
channel p^ ^ — » n^ 0,zi ^Tr^ is accessible, reducing 
the branching fraction of p^°^ — * Witt- 

At this point our data sample is still dominated by 
background. We therefore use additional variables that 
characterize the topology of the events to discriminate 
between signal and background. These variables are the 
azimuthal angle difference between the two jets A<p(j,j), 
the azimuthal angle difference between the electron and 
the fa, A<f>(e,$T), the transverse momentum of the dijet 
system pr{jj)j the scalar sum of the transverse momenta 



of the electron and the two jets Hj,, the invariant mass 
of the dijet system M(jj), and the invariant mass of the 
W boson-dijet system M(Wjj). The technicolor parti- 
cles are expected to have narrow widths (~ 1 GeV). We 
should therefore see enhancements in the distributions of 
M(jj) and M(Wjj), consistent in width with the detec- 
tor resolution. M(jj) corresponds to the reconstructed 
7Tt mass and M(Wjj) corresponds to the reconstructed 
Pt mass. We reconstruct the W boson from the electron 
and the missing transverse momentum using the W bo- 
son mass constraint to solve for p z of the neutrino. If 
there are two real solutions, we take the smaller value of 
neutrino \p z \. If there is only a complex solution, we take 
the real part. Distributions of these variables are shown 
in Figs. [1] and [2j We use two approaches to separate sig- 
nal and background, a cut-based analysis and a neural 
network analysis. 

The cut-based analysis is optimized using Monte Carlo 
simulations to maximize the ratio S/ \J~B for every set of 
technimeson mass values. S is the expected number of 
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FIG. 1: Distributions of A<f>(jj), pr(jj) and H T after final kinematic selection. The Witt signal is shown for M(pr) = 210 GeV 
and M(jit) = 110 GeV. Arrows at the bottom indicate the cuts applied in the cut-based analysis for the signal mass point 
shown. 
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FIG. 2: Distributions of M(Jj) and M(Wjj) after final kine- 
matic selection. The Witt signal is shown for M(pr) = 
210 GeV and M(ir T ) = 110 GeV. Arrows at the bottom indi- 
cate the cuts applied in the cut-based analysis for the signal 
mass point shown. 



Wttt events and B is the expected number of background 
events. For each topological variable, the S/y/B ratio is 
evaluated as a function of the value of the variable to 
determine a set of lower, upper, or window cuts which 
maximizes this ratio. 

The neural network analysis uses the topological vari- 
ables Hj,, A<f>(e,fir), A<fi(jj), pr(ii), the transverse mo- 
menta of both jets and of the electron and j>T- A two- 
stage neural network based on the Multi Layer Percep- 
tron algorithm [13] is used. The first stage consists 
of three independent networks which are trained to re- 
ject the three main backgrounds, top quark production, 
W + bb production, and all other W+jets production in- 
cluding heavy flavors. Each of these three networks has 
eight input nodes and one hidden layer with 24 nodes. 
The second stage network has three input nodes, con- 
nected to the outputs of the three networks in the first 
stage, and one hidden layer with six nodes. The sec- 
ond stage network is trained using all nine physics back- 



ground processes. The networks are trained separately 
for each set of technicolor mass values. We then apply 
the trained neural networks to the collider data, techni- 
color signals, and physics and instrumental backgrounds 
to obtain the discriminator output spectra. We opti- 
mize the discriminator cut for every set of techniparticle 
masses to maximize S/ \T~B~. 

There is no excess in our data over the expected back- 
ground. We compute upper limits on the pr — > Witt — > 
ev 66(c) production cross section times branching frac- 
tion. In the cut-based analysis, which is a simple count- 
ing experiment, we compute an upper 95% C.L. limit 
on the signal using Bayesian statistics [23]. The neural 
network analysis performs a maximum likelihood fit of 
the data in the M(pr), M(ttt) plane to signal and back- 
ground expectations. The backgrounds are constrained 
to their expected values within statistical and systematic 
uncertainties. The uncertainties in the background event 
yields total to 10-12% and the uncertainty in the sig- 
nal selection efficiency is 10% for the cut-based analysis 
and 20% for the neural-net based analysis. The largest 
contributions to the systematic uncertainties are due to 
jet reconstruction efficiency, jet energy scale, 6-tagging 
efficiency, and, only for the signal, from the difference 
between fast and fully simulated detector Monte Carlo. 
The 95% C.L. upper limit on the signal cross section is 
then determined by the number of signal events below 
which lies 95% of the integral over the resulting likeli- 
hood function. 

The expected sensitivity and the regions excluded at 
95% C.L. by both analyses in the M{pr), M(ttt) plane 
for M v = 500 GeV are shown in Fig. GO For M v 
100 GeV, only a small region around M{px) = 190 GeV 
and M(ttt) = 95 GeV can be excluded. We note from 
Fig. [3]^ a), that the expected sensitivity of the neural net- 
work analysis is better than that of the cut-based analy- 
sis, as indicated by the larger 95% C.L. exclusion region. 
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FIG. 3: Expected region of exclusion (a) and excluded region 
(b) at the 95% C.L. in the M(p T ), M{-k t ) plane for p T -> 
Witt — * ev bb(c) production with Mv = 500 GeV. Kinematic 
thresholds from Wnr and -kt^t are shown on the figures. 



We quote the observed 95% C.L. exclusion region in the 
M(p T ), M(tt t ) plane in Fig. Mb) by the neural network 
analysis as our measurement [2l|. 

The results presented in this paper cannot be com- 
pared directly to those previously published [23}. The 
CDF experiment did not use Ref. Q and 0, but rather 
the models described in the earlier paper of Ref. 
a precursor to the TCSM. The LEP experiments used 
Ref. Q in which the cross sections, while appropriate 
for narrow px production in qq collisions, are incorrect 
for off-resonance production in e + e~ collisions such as at 
LEP (see Ref. [24{). Although differences in the employed 
TC models preclude a direct comparison with previous 
searches, the current search achieves a higher sensitivity 
to the considered physics process. 
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